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ABSTRACT 

The existence of a classical bulge in disk galaxies holds important clue to the assembly history 
of galaxies. Finding observational evidence of very low mass classical bulges particularly in barred 
galaxies including our Milky Way, is a challenging task as the bar driven secular evolution might bring 
significant dynamical change to these bulges alongside the stellar disk. 

Using high-resolution N-body simulation, we show that if a cool stellar disk is assembled around 
a non-rotating low-mass classical bulge, the disk rapidly grows a strong bar within a few rotation 
time scales. Later, the bar driven secular process transform the initial classical bulge into a flattened 
rotating stellar system whose central part also have grown a bar-like component rotating in sync with 
the disk bar. During this time, a boxy/peanut (hereafter, B/P) bulge is formed via the buckling 
instability of the disk bar and the vertical extent of this B/P bulge being slightly higher than that of 
the classical bulge, it encompasses the whole classical bulge. The resulting composite bulge appears to 
be both photometrically and kinematically identical to a B/P bulge without any obvious signature of 
the classical component. Our analysis suggest that many barred galaxies in the local universe might 
be hiding such low-mass classical bulges. We suggest that stellar population and chemodynamical 
analysis might be required in establishing the evidence for such low-mass classical bulges. 

Subject headings: galaxies: bulges - galaxies:kinematics and dynamics - galaxies: structure - 
galaxies:evolution - galaxies:spiral, galaxies:halos 


1. INTRODUCTION 

Classical bulges are typically thought to have formed 
as a result of violent mergers fe auffnia nn et al.l 1 199,11 
iBaugh et al.111996I iHookins et al.l ^0911 or collapse of 
primordial gas clouds ( Eggen et al.lll96A or coalescence 
of giant clumps in high-red shift galaxies (llmmeli et al.l 
|200^jElniegreene^_aL l2008li or multiple mi nor mergers 
( Bournaud et al.l 12007 : iHopkins et al.l I2010H . Size and 
mass of these classical bulges depend on the process that 
formed these stellar system but perhaps little on the sub¬ 
sequent evolution. Observation suggests that the classi¬ 
cal bulges are generally dispersion dominated, spheroidal 
systems which might remain dynamically unchanged over 
several billion years. A number of recent studies suggest 
coexistence of classical bulges an d bars in disk galaxies 
(jGadottill2009l lErwin et al.ll201,'T[l . Of particular interest 
are the low-mass classical bulges in barred galaxies, such 
as our own Galaxy. N-body modelling of BRAVA stellar 
kinematics reveals that our Milky Way mig ht harbor a 
low-ni a ss classical bulge i nside it’s B/P bulge IShen et al.l 
(j20ir)l l: [Pi Matteo et al.l (j2015[) . For a current summary 
of the Milky W ay’s bulge, the reade rs are referred to 
iGerhardI (j2014l l. Recent simulation bv iSaha et al.l (j2012f ) 
showed that such a low-mass classical bulge is signifi¬ 
cantly modified during the secular evolution along with 
the cool stellar disk. It remains unclear whether such 
low-mass classical bulges can reliably be detected in ob¬ 
servation. 

In the initial phase of galaxy evolution, disk galax- 
ies (at high redsh ift) go through violent instabilities 
(|Genzel et all 1200611 and possibly have formed various 
non-axisymmetric structures (most prominent of which 
are bar and spiral arms) which, later, drive the slow 


secular evolution of these galaxies. In the secular 
phase, the most efficient way a disk galaxy evolves 
is through forming a bar which facilitates the redis¬ 
tribution of energy and angular momentum between 
the disk, d ark matter halo and the preexisting clas¬ 
sical bulge (iD ebattist a fc SellwoodI I2000I lAthanassoulal 
l2f)0il:ISaha et aLll2012ll . As the bar becomes stronger, it 
goes through buckling instability and form boxy/peanut 
bulges a s demonstrated in numerous N-body simulation 


studies (IGombes & SandersI Il981t IPfenniger & Norman 

199C 

: IRaha et al.lll991l 

Martinez-Valmiesta & Shlosman 

2004 

: ISaha et al.l 12012 

). More than 40% of all ob- 


served bulges are B/P bulges (|Lutticke et al.l I2000H 
and their observed p roper ties are nicely summarized in 
iLaurikainen fc Salol (l2015f) . The formation of a B/P 
bulge in the central region is an energetically favor¬ 
able phenomenon, as the B/P bulge is comparatively 
more random motion dominated than the initial disk. 
In essence, secular evolution transforms the central part 
of a cool stellar disk to one w ith a comparatively hot 
centr ally concentrated bulge (iKormendv fc KennicuttI 
1200411 . When such a process takes place in the presence of 
an initially non-rotating low-mass classical bulge at the 
center, the secular process also transforms the low-mass 
classical bulge to a fast rotating object whose kinematic 
properties resemble that of a B/P bul ge, namely havin g 
cylindrical rotation in the inner region (|Saha et al.ll2012fl . 
The final composite bulge is a superposition of the B/P 
bulge and the transformed classical bulge, both of which 
are cylindrically rotating. The goal of this paper is to 
find out whether such a low-mass classical bulge can be 
unmasked from the composite bulge. We present a sys¬ 
tematic morphological and kinematic analysis to extract 
information about the presence of the classical bulge in 







































































2 


Saha 


the barred galaxy. 

2. GALAXY MODEL AND SIMULATION 

The initial galaxy model, in equilibrium, consists of 
an initially axisymmetric exponential disk, a cored dark 
matter halo and a non-rotating classical bulge. The ini¬ 
tial classical bulge (hereafter, ICB) is strongly flattened 
by the disk gravity. The ICB in this galaxy model has 
a total mass of Mb = 0.06Md, Md being the disk mass, 
and initial ellipticity in edge-on projection is given by 
Cb = 0.46. The initial stellar disk is cool with Toomre 
Q = 1.4 at i? = 2.5Rd. F urther details on the model is 
given in iSaha et al.l (|2012ll . 

We scale the model such that Md = 4.58 x 
and Rd = 4.0 kpc. Then the time unit is given by 24.9 
Myr. We have used a total of 10 million particles. The 
softening lengths for the disk, bulge and halo particles 
used are 12, 5 and 20 pc respectively following the sug¬ 
gestion of lMcMillan fc DehnenI (l2007h. The simulation is 
performed using the Gadget code (iSoringel et all 1200 111 
with a tolerance parameter 9toi = 0.7, integration time 
step ~ 0.4 Myr. The simulation was evolved for a time 
period of ^ 3.0 Gyr. 

3. MORPHOLOGY OF THE COMPOSITE BULGE 

The initial axisymmetric stellar disk, being cool, 
rapidly forms a bar with two-armed spiral as shown in 
Fig. [J The spiral produces rad ial heating (jSaha et al.l 
l20ir)n^llwood fc Carlberd 120141 1. which eventually de¬ 
stroys itself leaving a strongly barred galaxy at the end 
of the simulation. The end product resembles a typ- 
ical barred lenticu lar galaxy without any spiral arm 
(|Gortesi et al.ll20l'^ . In the center of this barred galaxy 
lies the low-mass classical bulge (see rounder contours at 
t = 0 in Fig. [T]). However, as time progresses, even at 
the very central region (< 0.5Rd) the density contours 
are no longer rounder but are purely elliptical suggesting 
non-existence of the classical bulge component. We have 
checked this further inside R < 0.2Rd and found similar 
result. In fact, the face-on surface density maps of the 
classical bulge stars alone (bottom panel of Fig. [T]) show 
that the inner part of this bulge is actually a bar - to¬ 
gether with the disk stars, the central part of the galaxy 
behaves like a single bar component without any struc¬ 
tural signature of the classical bulge. It has been shown 
previously that such a classical bulge actually absorbs a 
significant fraction of the disk angular momenta emit¬ 
ted by the bar and is transformed into a rapidly rotating 
bar-like object (|Saha et al.ll2012D . 

In edge-on projection, as shown in Fig. [31 the galaxy 
model at t = 0 has a classical bulge which extends upto 
about 0.5Rd in the vertical direction along the minor 
axis. As the bar undergoes buckling instability, it forms 
a B/P bulge in the inner part of the disk. The outcome 
of this physical process is a composite bulge which is a 
superposition of the preexisting classical bulge and the 
B/P bulge. Such composite bulges are not uncommon 
and they are repor ted in a numb e r of barred ga.laxies , 
recently studied bv iGadottl (|2009l l: lErwin et al.l (j20I5l l. 
The vertical extent of the B/P bulge is about O.TRd (see 
Fig. |3|) which completely encompass the flattened classi¬ 
cal bulge. In other words, the classical bulge camouflages 
inside the B/P bulge. This may surprise some readers; a 
priori it is not clear what would be the vertical extent of a 


B/P bulge with a given size and mass of a classical bulge 
to begin with the initial disk. In any case, the generality 
remains unclear at the moment, as more such N-body 
experiments are needed. Obviously, massive and bigger 
classical bulges can be excluded as those can not remain 
immersed in the B/P bulge. Note the surface density 
maps of these composite bulges, shown in Fig. [31 appear 
identical to a pure B/P bulge with no morphological sig¬ 
nature of the low-mass classical bulge. In other words, 
the present study suggests that a B/P bulge in a galaxy 
might hide a low-mass classical bulge which carries valu¬ 
able information on the galaxy formation history. 

To confirm, we analyzed the surface density profiles 
at different epochs during the secular evolution of this 
galaxy model. Fig. [2] depicts the density profile taken at 
t = 113. We have performed 3-component decomposi¬ 
tion on this barred galaxy model: an ex ponentia l stella r 
disk, sersic bulge and a sersic bar model lCTadottil ([200911 . 
Our three component decomposition shows that the cen¬ 
tral region can be well modelled by a pure sersic bulge 
with sersic index n = 1.8 without a need for a classi¬ 
cal bulge, atleast in photometric sense. The radius at 
which the bulge and disk densities are equal is given by 
Rbd = 0-83Rd- It appears that the final barred galaxy 
contains no observable morphological evidence of a clas¬ 
sical bulge in the bulge dominated region {R < Rbd)- 
In the next section, we analyze the kinematics of this 
composite bulge. 


4. KINEMATIGS OF THE COMPOSITE BULGE 

Classical bulges are kinematically hotter with rota¬ 
tional motion intermediate between the ellip ticals and 
B/P bulges ([Kormendv &: IllingworthI 1198211 . As per 
kinematics concerned, B/P bulges are further distin¬ 
guished from the classical b ulges, in that B /P bulges 
possess cylindr i cal rotation (iKormendy 119^ . but see 
[Williams et al.l ((201 Ill : lSaha fc Gerhardl ( 20I3ll for excep¬ 
tional cases. The kinematics of B/P structures includ¬ 
ing various projection effects in simulation of dis k galax ¬ 
ies are discussed in depth by Debattista et al.l (|2005ll : 
llannnzzi fc Athanassonlal ([2015 1. 

Fig. |3| depicts the 2D line-of-sight kinematic map (in 
edge-on projection) e.g., mean velocity and velocity dis¬ 
persion, before and after the formation of the B/P bulge 
in our model galaxy. While creating these maps, we have 
rotated the bar so that its major axis is aligned with 
the x-axis and restricted to one projection only. Clearly, 
the ICB is non-rotating. Once the B/P bulge is formed 
and evolved, the final composite bulge (classical -|- B/P 
bulge) shows clear cylindrical rotation. In addition to 
the 2D kinematic map, we check this using long-slit kine¬ 
matics at three different heights above the disk mid-plane 
{z = 0) taken at three different epochs during the secular 
evolution, see Fig. H) It is evident from the figure that 
the classical bulge has a P/cr = 0 at t = 0. The disk 
rotation velocity rises outward and the velocity disper¬ 
sion remains nearly flat in the inner region, as a result 
the local Vja rises linearly outward. Note that the lo¬ 
cal Vja for the slit-I {z/Rd = 0) goes to zero at i? = 0 
and this basically probes the differentially rotating stel¬ 
lar disk. At subsequent times, we have calculated the 
local Vja at Rbd = 0.83i?d and in all cases, we found 
Vja > 1.0 suggesting a kinematically cool bulge. At 
t = 54, the composite bulge rotates cylindrically. The de- 
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Fig. 1.— Upper panel: face-on stellar surface density (including both disk stars as well as classical bulge stars) maps in model RCG004 
at t=0, 54 and 113 (internal unit). Bottom panel: face-on surface density maps of the classical bulge stars alone at same epochs. 


gree of cylindr ical rotation becomes s tronger with times 
(as shown bv iSaha fc Gerhard l2013fl . see the panel at 
t = 113 of Fig. m The major axis velocity dispersion is 
nearly flat within the extent of the B/P bulge or even 
within Rm = 0.83i?d and exhibit a shallow decline be¬ 
yond that. The local V/a at Rtd for the composite bulge 
at t = 113 is close to 1.1 whic h can be considere d as a 
kinematically cool component (jErwin et al.ll20l'^ . But 
even otherwise, the kinematic signatures are definitely 
not of a classical component. So far all the kinematic 
diagnostics suggest that this composite bulge is actually 
a B/P bulge both morphologically and kinematically. It 
remains puzzling how to separate classical bulge stars 
that are hidden inside the B/P bulge. 

We then examined the velocity histograms at different 
locations within the B/P bulge region, more specifically 
within Rbd = O-SSRd- First, we looked at the velocity 
histograms along the minor axis of the bulge as shown in 



Fig. 2.— Radial surface density profile of the barred galaxy at 
t = 113 and its decomposition into a bulge+bar+disk. The sersic 
index of the bulge is n = 1.82 and the bar is modelled with a sersic 
index n = 0.35. 


Fig. [5] at three different epochs during the secular evo¬ 
lution. Each figure (in both the panels) consists of 3 
histograms : for classical bulge stars alone, disk stars 
alone and the composite bulge (classical -b B/P bulge) 
stars. Note that such an exercise is possible in our N- 
body simulations as we have a unique ID attached to each 
particle. As often the case, the bulge, whether classical, 
B/P or composite bulge, do not show any net rotation 
about the mi nor axis. As time progresses, the bar heat 
up the stars (|Saha et al.ll201C)l) as a result of which all 
the histograms fatten. Interestingly, at t = 113, all the 
histograms (due to the classical bulge stars, B/P bulge 
stars and of the composite bulge) are nearly identical; 
in other words the minor axis histograms bear no dis¬ 
tinguishable sign of the classical bulge. Next we looked 
at a small area {x/Rd = [0.5, 0.6], ly/i^dl = 0.05) along 
the major axis of the composite bulge - this is depicted 
in the lower panel of the figure. Initially, the classical 
bulge stars show no net rotation (histograms in red) and 
the disk stars show net rotation as per construction of 
the initial model. At t = 54 and II3, the classical bulge 
stars show net rotation (see non-zero mean of the veloc¬ 
ity histograms) and this nearly coincides with the B/P 
bulge as well as the composite bulge. Although there is 
a separate rotating stellar component (here, the classi¬ 
cal bulge) sitting inside the B/P bulge, the net velocity 
histograms of the composite bulge is almost identical to 
that of the B/P bulge and no detectable signatures of 
the rotating classical bulge. 

At this point, it is useful to recall that the density 
structures of the classical bulge shows that there is a bar 
(see Fig. |T] in the central region. T his was referred t o 
as the classical bulge-bar, reported in ISaha et al.l (|20I2D . 
This classical bulge-bar rotates in sync with the disk bar 
as can be inferred from position angle of the two of bars 
at different times during the evolution. As a result of 
this, there is a fine blend of the classical bulge stars and 
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Fig. 3.— 2D line-of-sight moment maps (edge-on projection with the major axis of the bar aligned parallel to x-axis) at three different 
epochs during the evolution. The upper panel denotes surface density, middle panel the mean velocity and the lower panel the velocity 
dispersion. 
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Fig. 4.— Long-slit stellar kinematics from the region covering the 
composite bulge: the upper panel shows the line-of-sight stellar 
velocity profiles at 3 different heights above the disk mid-plane 
during the evolution. The lower panel shows the same but for 
the velocity dispersion profiles. Note that at t = 0, the slit at 
zIRd = 0.1 and 0.22 essentially probe the classical bulge stars. 

disk stars that are part of the bar/B/P bulge, explain¬ 
ing the near identicalness of the velocity histograms. In 
other words, this would explain why the central region, 
although a composite bulge, behaves both morphologi¬ 
cally and kinematically like a B/P bulge with a perfectly 
hidden low-mass classical bulge. 
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Fig. 5.— Velocity histogram along the minor (upper panels) and 
major axis (lower panels) of the galaxy model. There is no partic¬ 
ular signature of the classical bulge component that resides inside 
the boxy/peanut bulge. Each histogram is fitted with a Gaussian 
distribution. Blue and red histograms are due to disk and classical 
bulge stars alone. The solid black lines represent the histograms of 
the composite bulge. 


5. DISCUSSION AND CONCLUSIONS 

It is generally believed that spheroidal stellar system, 
of which classical bulges are a subset, do not take take 
part in the secular evolution. Indeed, this may be true 
if such classical b ulges are massive and bigger. But it 
has been shown bv ISaha et al.l (I20I2D that very low mass 
classical bulges get substantially affected by the secular 
evolution driven by a bar. In fact, what bar driven sec¬ 
ular evolution does to a cool stellar disk is to produc e 
a central concentration (|Kormendv fc Kennicutdl2004D . 
essentially growing a hotter bulge with a V/a less than 
that of the disk. Kinematically, a B/P bulge is essentially 
a vertically thickened bar with hotter stars than what lies 
beyond the bar’s corotation. What secular evolution did 
to this low mass classical bulge is substantial and quite 
similar to th e disk. The centra l part of this classical bulge 
grew a bar (|Saha et al.ll20I2l) that is nearly identical to 


the disk bar. The bar later heat the bulge stars as is evi¬ 
dent from the minor axis velocity histograms. From this 
point of view, the classical bulge went through a secular 
evolution parallel to the disk with similar end products, 
although not quite the same. The generality of this out¬ 
come is not clear at the moment as it requires examining 
a large number of such high resolution N-body simula¬ 
tions and will be the subject matter of a future paper 
in the series. In any case, the current analysis produces 
a twist in a sense that an apparently pure B/P bulge 
could actually be a composite bulge with photometrically 
and kinematically hidden low-mass classical bulge which 
coul d possibly have their o rigin in mergers/minor merg¬ 
ers (|Khochfar fc Sillj I2006D in accordance with ACDM 
paradigm of galaxy formation. 

The next obvious question is whether the presence of 
these very low mass classical bulges can ever be estab- 















































Lost in secular evolution 


5 


lished unambiguously. Observationally, classical bulges 
tend to show older stellar population, quiescent in star 
formation activity and populate, in general, the red- 
sequence of the color-magnitude diagram. In contrast, 
pseudobulges/B/P bulges exhibit intense star forma¬ 
tion and tend to occupy the blue cloud of the color- 
magnitude diagram. Interestingly, galaxies with compar¬ 
atively lower mass classical bulges show some amount of 
star format ion activity an d about half of these galaxies 
are barred (|Gadottil 12009( 1. In other words, these barred 
galaxies would be hosting composite bulges with star for¬ 
mation activity intermediate between the pseudobulges 
and the ones hosting massive classical b ulges. If these 
star formation is drive n by minor mergers (|Kaviraill2014l : 
iSachdeva et al.1 (201511 . the preexisting classical bulge 
would further be contaminated and might increase the 
level of complexity for the stellar population analysis 
to draw a firm conclusion on the evidence of a pris¬ 
tine classical bulge. Metallicity gradient could be an im¬ 
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portant diagnostic to argue for the evidence of a clas- 
sical bulge, as it has been done for our Milky Way 
((Zoccali et al.[[2008l: (Gonzalez et al.l[2M^ (Johnson et al.l 

(201.111 . Even this can turn out to be a difficult task as the 
bar would act as a strong a gent for mixing and migration 
(jSellwood fc Binnevl (2002ll which would affect the clas¬ 
sical bulge stars, might as well erase its imprint. How¬ 
ever, if these bulges are formed at the very early epoch of 
the galaxy assembly history, they might contain extreme 
metal poor stars like those found in the Mhky W ay’s 
bulge (|Garcia Perez et id] (201.'K : (Howes et al.( (20141 and 
might convey an important message about the presence 
of an old classical bulge; of course a number of possibil¬ 
ity remains as they could be the halo stars or thick disk 
stars. In the end, one perhaps has to rely on chemody- 
namical analysis to obtain a reliable answer. 
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